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ABSTRACT 


A  FORTRAN  program  has  been  developed  to  allow  for  the 
use  of  probabilistic  design  methods  in  the  numerical  optimi¬ 
zation  process .  The  program  was  written  as  a  set  of  sub¬ 
routines  for  COPES  (Control  Program  For  Engineering  Synthesis) . 
COPES  maximizes  or  minimizes  a  numerically  defined  objective 
function  subject  to  a  set  of  inequality  constraints  using 
the  optimization  program  CONMIN  (A  Fortran  Program  for  Con¬ 
strained  Function  Minimization) .  The  program  developed  here 
allows  for  the  use  of  both  the  normal  and  lognormal  distribu¬ 
tion  models.  Design  examples  are  presented  to  demonstrate 
the  program  capabilities.  User  instructions  are  provided  for 
inclusion  in  the  COPES  user's  manual. 
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I.  INTRODUCTION 


In  the  traditional  deterministic  approach  in  engineering 
design  the  variables  are  treated  as  single  valued  numbers. 
These  variables,  whether  dimensions,  material  properties, 
loads,  etc. ,  are  actually  statistical  in  nature.  By  the  use 
of  safety  factors  designers  protect  against  these  variations 
by  usually  over  designing.  A  more  logical  approach  is  to 
take  into  account  the  statistical  data  known  about  each 
variable  and  design  for  a  certain  reliability.  By  consider¬ 
ing  the  statistical  nature  of  each  variable  we  should  be 
able  to  better  predict  reliability  and  performance. 

Numerical  optimization  has  proven  to  be  a  very  powerful 
tool  in  engineering  design.  Virtually  all  design  problems 
require  minimization  or  maximization  of  some  objective.  For 
the  design  to  be  acceptable,  it  must  also  satisfy  a  certain 
set  of  specified  requirements  called  constraints.  If  these 
constraints  are  specified  as  probabilities  of  failure  the 
result  would  be  a  probabilistic  design. 

The  purpose  of  this  research  was  to  test  the  applica¬ 
bility  of  combining  probabilistic  design  concepts  with  those 
of  numerical  optimization  by:  (1)  developing  a  pilot  computer 
code  to  calculate  probabilities  of  failure,  (2)  incorporate 
this  computer  code  with  COPES/CONMIN  in  order  to  perform 
probabilistic  optimization,  and  (3)  test  using  numerous 
examples . 


It  is  assumed  throughout  this  discussion  that  the  reader 
is  familiar  with  the  use  of  the  COPES/CONMIN  optimization 
program  [Ref.  1],  and  [Ref.  2]. 
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II.  NUMERICAL  OPTIMIZATION 

The  numerical  optimization  problem  considered  here  is 
stated  as  follows:  Find  the  set  of  n  design  variables 
contained  in  the  vector  X  which  will 

Minimize  F(X)  (1) 

Subject  to: 

gj  (X)  <_  0  j  =  l,m  (2) 

x^  <  x.  <  xV  i  =  l,n  (3) 

l-i-i 

The  components,  x^,  of  X  are  referred  to  as  design  variables 
which  are  changed  to  improve  the  design.  The  function  F(X) 
is  called  the  objective.  Inequality  constraints,  g^ (X) ,  are 
the  response  limits  imposed  on  the  design.  There  are  n  de¬ 
sign  variables  and  m  inequality  constraints.  The  lower  and 
upper  bounds,  x^  and  x^  are  limits  imposed  on  the  design 
variables  to  insure  a  practical  result. 

For  a  deterministic  design  COPES  [Ref.  1],  determines 
the  constraints  and  objective  as  follows.  Consider  the 
design  of  a  single  bar  undergoing  uniaxial  tension, 
wish  to  minimize  the  weight,  the  objective, 


If  we 


where 


p  =  the  specific  weight; 

A  =  the  area  of  the  bar; 

L  =  the  length  of  the  bar. 

The  design  variable  is  the  area,  A.  If  the  stress  in  the 
bar  is  S  and  the  stress  limit  imposed  is  Sy  we  desire  that 
S  be  less  than  or  equal  to  Sy  so  the  constraint  becomes 

S/Sy  -  1  <  0  (5) 

In  this  investigation  a  computer  subroutine  was  developed 
to  provide  constraints  based  on  probabilities  of  failure. 

If  the  allowed  probability  of  failure  is  PF  and  the 
actual  probability  of  failure  was  calculated  to  be  Pf,  then 
the  constraint  would  be 

Pf/PF  -  1  <  0  (6) 


III.  PROBABILISTIC  DESIGN 


A.  RANDOM  VARIABLES 

If  one  has  a  large  population  of  n  elements  for  which 
some  parameter  x^,i*l,n  is  determined,  the  mean  value  of 
the  population  is 


i  n 

i  -  -  I  (x.) 

x  n  j£l  1 


The  variance  is 


i  n 

V  =  -  l  (x.  -  y  ) 

x  n  1  Mx 


(7) 


(8) 


The  standard  deviation  is 


o 


X 


/v- 

X 


(9) 


The  standard  deviation  is  a  measure  of  the  amount  of  varia¬ 
bility  of  the  population  data.  An  additional  measure  of 
the  variability  is  the  dimensionless  coefficient  of  variation. 


C 


x 


(10) 


In  engineering  design  the  coefficient  of  variation  is  fre¬ 
quently  used. 
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B.  FUNCTIONS  OF  SEVERAL  VARIABLES 


If  Y  is  a  function  of  n  variables,  the  standard  deviation 
of  the  function  Y  can  be  approximated  as  follows: 

If  Y  =  ft*.;  i  *  l,n) 

Assuming  that  all  are  independent  random  variables 
then  according  to  [Ref.  3,  p.  59] 


[  l  °l  11/2 

i=l  dXi 


x . 
1 


(11) 


C.  STATISTICAL  MODELS 

1 .  Normal  Distribution 

The  normal  distribution  model  is  one  of  the  most 
widely  used.  The  normal  density  function  is  given  by 


where : 


fx(x) 


1  ■(x"wx) 

- _  exp[ - ^ - ] 

/2tto  2  a 


(12) 


fxU) 


the  probability  density  function  of 
the  random  variable  X; 


standard  deviation  of  the  random 
variable  X; 


»  mean  of  the  random  variable  X. 
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The  magnitude  of  the  standard  deviation  determines  the  dis¬ 
persion  of  the  distribution.  A  large  standard  deviation  will 
result  in  a  wide  bell  curve  and  a  small  standard  deviation 
will  result  in  a  narrow  one.  The  curve  is  symmetric  about 
the  mean.  Approximately  68.26%  of  all  samples  will  lie 
within  plus  or  minus  one  standard  deviation  of  the  mean  value. 
The  standardized  form  of  the  normal  density  function  [Ref.  4: 
p.  194],  can  be  expressed  as 


$Y(y) 


7F  exp  t=j£i 

✓  2it 


(13) 


where 


y  = 


(14) 


<t>Y(y)  =  the  standardized  form  of  the  single 

variate  normal  density  function. 

The  single  variate  cumulative  normal  distribution  function 
is  given  by 


X 

Fx(x)  =  /  fx(x)dx 


(15) 


where : 


Fx(x)  *  the  probability  that  the  random  variable 
*  X  is  equal  to  or  less  than  the  specific 

value  x. 
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Or,  using  the  expression. 


x  -  u. 


y  = 


Equation  (15)  can  be  rewritten  as,  [Ref.  4:  p.  195] 


PX(X) 


/2ir 


<x'*'x)/ox 


exp [— *— ] dy  (16) 


or 


F„(x)  =  $v  [  ■  ]  (17) 

X  Y  °x 

2 .  Lognormal  Distribution 

If  ln(X)  has  a  normal  distribution  then  X  is  said  to 
have  a  lognormal  distribution.  The  lognormal  density  func¬ 
tion,  [Ref.  4:  p.  196],  is  given  by 


£x(x) 


-^|exPr(ln(*1f> 

/2’°x  2ox 


-] ;  x  >  0 


(18) 


where : 

A 

X  =*  mean  of  ln(X)  ; 

A 

ax  =  standard  deviation  of  ln(X) . 


The  cumulative  lognormal  distribution  function,  [Ref.  4: 
p.  196],  is  given  by 

ln(x) -X 

A 

^  x  2 

Fy (x)  «  —  /  exp[3f-]dy;  x  >  0  (19) 


or 


Fjt(X)  .  *Yt 


(20) 


The  mean,  yx,  and  the  variance  v  ,  of  X  are 


a2 

x, 


y  x  ■  exp[X  +  ~y] 


(21) 


Vx  *  exp[2X + a2] [exp(a2)  -  1] 


(22) 


Solving  for  X  and  a. 


X  =  In (y 2)  -  |  In  (Vx  +  y2) 


(23) 


^/-ln(y2)  +  ln[Vx+y2] 


(24) 


D.  PROBABILITY  OF  FAILURE 

Consider  the  constant  area  bar  in  Figure  3.1  under  load 
P.  If  the  stress  is  the  bar  is  S  and  the  stress  limit  in 
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Figure  3.1.  BAR  UNDER  UNIAXIAL  LOAD 

tension  for  the  material  is  Sy,  then  the  probability  of 
failure  of  the  bar  is  the  probability  that  the  stress  S  will 
be  greater  than  the  limit  Sy.  The  method  for  determining 
the  probability  of  failure  in  the  bar  for  both  the  normal 
and  the  lognormal  formats  is  presented  below. 

1 .  Normal  Distribution 

If  both  S  and  Sy  are  random  variables  which  follow  a 
normal  distribution  then  the  probability  of  failure  is  the 
probability  the  S  is  greater  th an  Sy.  The  failure  function, 
[Ref.  5:  p.  8-24] 


Z  =  (Sy  -  S) 


(25) 
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is  also  normally  distributed.  The  probability  of  failure 
is  the  probability  that  the  failure  function  is  less  than 
or  equal  to  0.0  or 


Pf  ■  P(Z<_0) 

In  the  standard  form 

Pf  =  virf1 

Where,  by  the  algebra  of  functions,  [Ref.  3] 


(26) 


(27) 


(28) 


and 


a 


z 


(29) 


The  safety  index,  [Ref.  5:  p.  8-24],  is  defined  as 


S 


The  probability  of  failure  is  then 


(30) 


pf  =  4gi-ei 


(31) 
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or  the  area  under  the  unit  normal  density  function  from  8 

tO  oo  . 


Pf  =  /  —  exp  [-4—]  dg  (32) 

8  /sr  2 

2 .  Lognormal  Distribution 

If  the  same  two  variables  S  and  Sy  follow  a  logri- 
thetic  normal  distribution  then  the  probability  of  failure 
is  expressed  by,  [Ref.  5:  p.  8-271 

Pf  »  P(Sy<S)  =  P(ln(Sy/S)  <  0)  (33) 


Letting 


Z  =  In (Sy/S) 


(34) 


then 


Pf  =  P(Z  <  0) 


(35) 


or 


Pf  =  *3t-81 

The  safety  index,  [Ref.  5:  p.  8-29]  is 


(36) 
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ln(Sy/S) 


JlMuZ 


CSy»(1  +  Csn 


where  the  Sy  and  S  are  median  values  and 


(37) 


where 


(38) 


(39) 


(40) 

(41) 


6  will  be  normally  distributed  and  the  probability  of 
failure  is 


Pf 


[ 


-ln(Sy/S) 


4 


In [ ( 1+Cgy) ( 1+Cg) ] 


] 


(42) 


The  probability  that  S  is  greater  than  Sy  is  then  determined 
in  the  same  manner  as  the  normally  distributed  case,  by 
integrating  the  unit  normal  density  function  from  S  to  ®. 
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IV.  PROBABILISTIC  OPTIMIZATION 


The  objective  of  probabilistic  optimization  is  to  mini¬ 
mize  or  maximize  the  mean  value  of  a  given  function  subject 
to  a  set  of  constraints  based  on  allowed  probabilities  of 
failure. 

A.  CONSTRAINTS 

The  computer  subroutine  developed  in  this  investigation, 
C0PE19,  calculates  the  constraints  as  described  in  Chapter 
II.  Given  a  constrained  variable  in  terms  of  a  desired 
probability  of  failure,  one  needs  only  to  calculate  the 
actual  probability  of  failure  to  determine  the  constraint 
value . 

Consider  the  cantilevered  beam  in  Figure  4.1.  A  proba¬ 
bility  of  failure,  PF,  is  desired  for  stress  in  bending. 

The  mean  values  and  coefficients  of  variation  for  the  yield 
strength,  Sy,  the  length  L,  the  load  P,  and  the  dimensions  B 
and  H  are  given.  The  stress  in  bending 


The  standard  deviation  of  the  bending  stress 


(43) 


3SB 


SB 


l'  '  3H 


,2.  3SB  T" 

JH  ^TeTV 


(44) 
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The  safety  index 


6 


The  calculated  probability  of  failure 

Pf  =  (46) 

The  constraint  will  be 

G  =  ||  -  1  (47) 

In  the  computer  program  developed,  C0PE19,  the  required 
partial  derivatives  of  the  constrained  variables  with 
respect  to  the  variables  in  the  expression  are  calculated 
by  finite  difference  methods. 

B.  OPTIMIZATION 

Consider  the  three  bar  truss  in  Figure  4.2.  The  geometry 
is  specified  and  there  are  two  independent  load  conditions. 
The  design  task  is  to  determine  the  areas  required  for  each 
bar  that  will  yield  a  minimum  mean  value  of  the  structure 
weight  while  keeping  the  stress  in  any  member  less  than  the 
yield  stress. 


Sy  -  SB 


(45) 
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THREE  BAR  TRUSS 

H - e+3 - 


ll\v 


Figure  4.2.  THREE  BAR  TRUSS 
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Minimize:  W  =  pH(/2Al  +  A2  +  /?A3) 

Subject  to:  Sy  compression  <_  SIG^j  £  Sy  tension 
where  SIGi ^  is  the  stress  in  member  i  under  load  condition  j . 
Given:  Material  p  =  .1  lb./cu.  in. 

Geometry  H  =  10 . 0  in . 

Loads  PI  =  P2  =  20,000  lb. 

Stress  limits:  Sy  compression  =  -15,000  psi 

Sy  tension  =  20,000  psi 

Maintain  symmetry  Al  =  A3 

Beginning  with  the  design  Al  *  A2  *  A3  =  1.0  sq.  in.  and 
optimizing  using  COPES/CONMIN  (deterministic)  the  following 
results  are  obtained: 

W  =  2.632  lb. 

Al  =  A3  =  .7796  sq.  in. 

A2  =  .4275  sq.  in. 


The  critical  constraints  are  SIGH  and  SIG32,  which  are 
equal.  The  design  space  is  shown  in  Figure  4.3.  As  no 
safety  factor  was  used,  the  above  deterministic  design 
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Figure  4.3.  THREE  BAR  TRUSS  DESIGN  SPACE  (DETERMINISTIC) 
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corresponds  to  a  probability  of  failure  of  0.5  or  a  50% 
chance  that  SIGH  and  SIG32  will  be  greater  them  the  stress 
limit  Sy.  Utilizing  the  subroutine  COPE19  to  provide  proba¬ 
bilistic  based  constraints  and  assuming  that  the  coefficients 
of  variation  of  all  the  variables  are  0.1,  the  following 
results  are  obtained  when  designing  for  a  maximum  probability 
of  failure  of  0.5  in  each  member?  assuming  a  normal  distribution, 

W  «  2.639  lb. 

Al  =  A3  =  .7883  sq.  in. 

A2  =  .4094  sq.  in. 

The  design  space  is  shown  in  Figure  4.4.  The  lines  of  con¬ 
stant  probability  of  failure  for  SIGll  are  shown.  As  one 
would  expect  the  line  of  constant  0.5  probability  of  failure 
for  SIGll  in  Figure  4.4  is  the  same  as  the  constant  20,000 
psi  line  for  SIGll  in  Figure  4.3.  The  reverse  is  also  true. 

If  the  three  bar  truss  is  optimized  for  a  maximum  probability 
of  failure  of  0.4  the  resulting  safety  factor  for  SIGll  is 
1.03  or  19,361  psi  and  the  minimum  weight  is  2.726  lb.  Using 
19,361  psi  as  the  limit  of  stress  for  a  deterministic  optimi¬ 
zation  yields  approxiamtely  the  same  minimum  weight.  In  both 
cases  the  feasible  design  space  is  the  same  with  the  line 
of  constant  0.4  probability  of  failure  and  constant  19,361 
psi  lying  on  top  of  each  other. 
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R1  (SQ.  IN.) 

—  CONSTANT  OBJECTIVE  (LBS) 

-  CONSTANT  PROBABILITY  OF 

FAILURE  SIGH 

Figure  4.4.  THREE  BAR  TRUSS  DESIGN  SPACE  (PROBABILISTIC) 
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Table  I  shows  the  results  of  minimizing  the  weight  of 
the  three  bar  truss  for  various  probabilities  of  failure. 
Table  II  provides  the  deterministic  optimum  designs  for  the 
equivalent  factor  of  safety.  The  difference  in  designs  is 
within  the  numerical  accuracy  of  the  optimization. 


TABLE  I 

THREE  BAR  TRUSS  PROBABILISTIC  DESIGNS 


PROBABILITY 

WEIGHT 

A1  =  A3 

A2 

SIGH 

SAFETY 

OF  FAILURE 

LBS. 

SQ.  IN. 

SQ.  IN. 

PSI 

FACTOR 

.5 

2.6390 

.78827 

.40938 

20000 

1.00 

.4 

2.7260 

.81606 

.41785 

19361 

1.03 

.3 

2.8238 

.84478 

.43438 

18691 

1.07 

.2 

2.94  30 

.86887 

.48544 

17940 

1.11 

.1 

3.1208 

.92482 

.50500 

16914 

1.18 

.01 

3.6099 

1.0490 

.64280 

14640 

1.37 

.001 

4.0459 

1.1777 

.71493 

13060 

1.53 

.0001 

4.4770 

1.3212 

.73994 

11792 

1.70 

TABLE 

II 

THREE  BAR 

TRUSS  DETERMINISTIC 

DESIGNS 

SAFETY 

WEIGHT 

Al  =  A3 

A2 

SIGll 

FACTOR 

LBS. 

SQ.  IN. 

SQ.  IN. 

PSI 

1.00 

2.6326 

.77962 

.42752 

20051 

1.03 

2.7212 

.80145 

.45437 

19402 

1.07 

2.8200 

.  82773 

.47889 

18726 

1.11 

2.9322 

.87045 

.47014 

18002 

1.18 

3.1180 

.91845 

.52000 

16933 

1.37 

3.5936 

1.0807 

.53693 

14688 

1.53 

4.0258 

1.2033 

.62227 

13110 

1.70 

4.4800 

1.3379 

.69573 

11781 

It  should  be  noted  that  the  correspondence  of  the  proba¬ 
bility  of  failure  and  the  safety  factor  shown  here  is  unique 
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to  this  problem,  in  which  all  C  -  0.1.  In  general  nonlinear 
optimization,  this  one  to  one  relationship  cannot  be  assured. 

C.  EFFECTS  OF  COEFFICIENTS  OF  VARIATION  ON  THE  DESIGN  SPACE 

As  the  coefficients  of  variation  of  the  variables  used  in 
the  design  subroutine  become  smaller,  the  bands  in  the  design 
space  for  the  constrained  values  between  Pf  =  1.0  and  Pf  =  0.0 
becomes  narrower.  Figure  4.5  shows  the  design  space  for  the 
three  bar  truss  in  the  vicinity  of  the  optimum  for  the  case 
where  all  variables  have  coefficients  of  variation,  C  =  0.1. 
Figure  4.6  shows  the  same  space  for  the  case  where  all  coeffi¬ 
cients  of  variation,  C  =  0.001.  In  both  figures  the  lines  of 
constant  objective  and  probabilities  of  failure  for  the 
critical  constrained  variable,  SIGH,  are  shown.  The  line 
of  constant  Pf  *  0.5  is  the  same  for  both  but  the  band  between 
Pf  -  0.9  and  Pf  =  0.1  is  much  narrower  for  the  case  where 
C  =  0.001. 

This  phenomenon  can  cause  numerical  difficulties  in  the 
optimization  process.  Premature  termination  of  optimization 
can  occur  as  the  optimizer  may  be  unable  to  move  down  the 
resulting  'narrow  valley'.  The  absolute  and  relative  termina¬ 
tion  criteria  of  CONMIN  may  have  been  satisfied  but  the  true 
optimum  might  not  have  been  reached. 

If  the  initial  design  is  a  considerable  distance  from 
the  narrow  band  of  changing  probabilities  of  failure,  CONMIN 
may  obtain  zero  gradients  for  the  constraints  which  would 
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Figure  4.5.  THREE  BAR  TRUSS  DESIGN  SPACE  (C  =  0.1) 
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Figure  4.6.  THREE  BAR  TRUSS  DESIGN  SPACE  (C  =  0.001) 


also  cause  termination  of  the  optimization  process.  A 
routine  to  attempt  to  prevent  this  was  included  in  C0PE19. 
The  cumulative  distribution  curve  for  the  normal  format  is 
shown  in  Figure  4.7.  The  numerical  integration  routine  used 
to  calculate  the  area  under  the  unit  normal  distribution 
curve  yields  a  probability  of  failure,  Pf  of  1.0000000  for 
a  8  of  -5.2999821  and  a  Pf  of  0.0000000  for  a  8  of  5.5000086. 
It  was  modified  to  "widen"  the  zone  between  Pf  =  1.0  and 
Pf  =  0.  A  polynomial  (a  line  equation)  was  placed  at  each 
end  of  the  cumulative  distribution  routine  to  yield  a  Pf 
of  1.0  at  3  equals  -25.0,  decreasing  to  a  Pf  of  0.9999866 
at  -4.2  and  a  Pf  of  0.00001335144  at  8  equals  4.2  decreasing 
to  a  Pf  of  0.0  at  25.0.  This  is  shown  in  Figure  4.8.  This 
provides  the  same  optimization  results  within  the  accuracy 
of  the  program,  but  avoids  the  problem  of  zero  gradients  of 
the  constrained  variables. 

An  additional  approximation  was  used  to  determine  the 
probability  of  failure  when  the  coefficients  of  variation 
concerned  are  zero  or  nearly  so.  The  probability  that  X  is 
greater  than  Y  is,  in  the  standard  format 

v-03  =  v- ^3  (48) 


where 


u 


z 


(49) 
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Figure  4.8.  MODIFIED  CUMULATIVE  NORMAL  DISTRIBUTION  CURVE 
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(50) 


and 


Mathematically,  if  a  =0,  then  one  of  three  situations  can 

z 

occur. 

1.  X  ■  Y,  6  =  0,  then  Pf  =  0.5; 

2.  X  >  Y ,  0  ,  then  Pf  =  1.0; 

3.  X  <  Y,  3  -  then  Pf  =  0.0. 

This  is  shown  graphically  in  Figure  4.9.  In  order  to  provide 
CONMIN  with  a  smooth  function  with  which  to  work,  a  polynomial 
was  used  to  connect  the  Pf  =  1.0  and  Pf  =  0.0  in  Figure  4.9. 
This  is  shown  in  Figure  4.10.  The  interpolating  polynomial 
is  used  when  J,  the  modified  safety  index,  is  within  one 
unit  of  0.  Where  3  is  defined  to  be 


3  = 


Y  -  X 


(51) 


A  similar  approximation  was  developed  for  the  lognormal 
model. 

Experience  has  shown  that  these  modifications  dramatically 
improve  the  numerical  stability  of  the  optimization  process 
without  significantly  affecting  the  mathematical  result. 
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V.  COMPUTER  PROGRAM 


A.  COPES  INPUT  DATA 

Standard  deterministic  optimization  requires  that  COPES 
data  blocks  A-Q  and  V  be  used.  For  probabilistic  optimiza¬ 
tion  the  same  data  is  entered  in  these  blocks  with  the 
exception  of  NCALC  and  IPROB  in  data  block  B.  Here  an  NCALC 
=  7  would  be  entered  to  indicate  probabilistic  optimization 
was  to  be  performed  and  an  IPROB  =  1  or  2  to  indicate  which 
probability  model  is  to  be  used.  Additional  data  blocks  are 
required  for  probabilistic  optimization.  Appendix  A  con¬ 
tains  detailed  input  instructions  for  these  data  blocks  as 
well  as  an  example  problem  including  sample  input  data. 

B.  OPTIMIZATION  DATA  FLOW 

The  interaction  between  COPES,  CONMIN,  and  the  analysis 
subroutine  will  be  the  same  for  probabilistic  optimization 
as  it  is  in  standard  deterministic  optimization  with  one 
exception.  In  the  optimization  process,  each  time  COPES 
requires  the  determination  of  a  constraint  vector  it  will 
call  COPE  19.  C0PE19  will  call  the  analysis  subroutine  as 
many  times  as  necessary  to  determine  the  required  partial 
derivatives  by  finite  difference  steps.  Once  these  are 
obtained,  the  constraint  vector  is  determined  and  provided 
to  COPES.  A  flow  diagram  for  probabilistic  optimization 
is  provided  in  Figure  5.1. 
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SUBROUTINE 
ANAL  1 2 


OUTPUT 


Figure  5.1.  PROBABILISTIC  OPTIMISATION  FLOW  DIAGRAM 


43 


C.  COPE 19  SUBROUTINE 

The  main  computer  subroutine  developed  in  this  investi¬ 
gation  is  C0PE19.  Three  additional  subroutines,  COPE20 
through  COPE22,  were  also  developed  to  perform  specific 
tasks  required  by  C0PE19.  The  operation  and  flow  of  these 
subroutines  are  discussed  below. 

1.  Input 

The  following  is  provided  to  C0PE19: 

a.  The  COPES  control  parameter  ICALC. 

b.  Which  probability  model  to  use,  either  normal 
or  lognormal. 

c.  The  number  of  probability  variables  which  make 
up  the  constraint  equations  (IVAR) . 

d.  The  global  locations  in  the  ANALIZ  subroutine 
common  block  of  each  IVAR. 

e.  Coefficient  of  variation  of  each  IVAR,  assumed 
constant  throughout  the  optimization. 

f.  The  number  of  constrained  variables. 

g.  The  global  location  in  the  ANALIZ  subroutine 
common  block  of  each  of  the  constrained 
variables . 

h.  The  upper  and  lower  limits  imposed  on  each  of 
the  constrained  variables. 

i.  The  allowed  probability  of  failure  at  the  upper 
and  lower  limits  of  each  constrained  variable. 


j.  The  coefficient  of  variation  of  the  constrained 
variable  limits. 


k.  The  variable  MGRAD  which  determines  how  often 

the  constrained  variable  gradients  are  calculated. 

2 .  Determination  of  Constraint  Vector 

When  COPES  requires  a  constraint  vector  for  a  given 
design  it  will  call  C0PE19  with  ICALC  =  2.  C0PE19  will  then 
calculate  the  partial  derivatives  of  each  constrained  varia¬ 
ble  with  respect  to  the  variables  which  make  up  the  con¬ 
straint  equations.  The  subroutine  COPE20  controls  the  finite 


difference  calculations 

• 

In  the  example  in 

Section  IV. A, 

the  cantilevered  beam, 

it 

would  calculate: 

3SB 

6P 

BH2 

(52) 

3L 

= 

3SB 

3P 

— 

6L 

BH2 

(53) 

3  SB 

12PL 

(54) 

3H 

= 

BH 

3SB 

6PL 

(55) 

3B 

2  2 

B  H 

C0PE19  will  then  calculate  the  standard  deviation  (Equation 
44)  of  the  constrained  variable,  in  this  case  SB,  as  dis¬ 
cussed  in  Section  III.B  and  Section  IV. A.  Now  the  safety- 
index  is  calculated  and  subroutine  C0PE21  called  to  deter¬ 
mine  the  probability  of  failure.  When  required,  C0PE21  will 
call  subroutine  COPE22  to  perform  numerical  integration  of 
the  normal  density  function.  The  constraint  value  to  be 
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stored  in  the  G  vector  is  then  determined  (Equation  6) . 

The  above  would  be  performed  for  each  constrained  variable 
and  the  resulting  G  vector  provided  to  COPES. 

D.  REDUCTION  OF  FUNCTION  EVALUATIONS 

In  order  to  determine  the  gradients  of  the  constrained 
variables  finite  difference  steps  are  taken  resulting  in  a 
significant  increase  in  the  number  of  analyses  performed  in 
the  optimization  process.  In  many  problems  these  gradients 
remain  essentially  the  same  throughout  the  optimization. 
Therefore,  it  is  sometimes  possible  to  reduce  the  frequency 
with  which  the  gradients  are  calculated  and  achieve  approxi¬ 
mately  the  same  result.  The  input  variable  code  LGRAD  is 
included  as  an  option  in  probabilistic  optimization  to  allow 
this  choice.  LGRAD  has  the  following  meanings: 

0:  Calculate  gradients  each  time  COPE19  is  called. 

1:  Calculate  gradients  at  the  beginning  of  each 
CONMIN  iteration. 

2:  Calculate  gradients  only  at  the  beginning  of 
optimization. 

Table  III  demonstrates  the  result  of  various  LGRAD  selections 
when  designing  the  three  bar  truss  for  an  allowed  probability 
of  failure  of  0.50,  where  all  coefficients  of  variations, 

C  *  0.1. 

Probab il is tic  optimization  was  performed  using  the  three 
LGRAD  values  on  numerous  additional  test  cases.  In  almost 
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TABLE  III 


REDUCTION  OF  FUNCTION  EVALUATIONS:  THREE  BAR  TRUSS  DESIGN 


LGRAD 

WEIGHT 

LBS. 

Al  =  A3 
SQ.  IN. 

A2 

SQ.  IN. 

NUMBER  OF 
FUNCT  EVALS 

Deterministic 

2.6326 

.77962 

.42752 

35 

Probabilistic 

0 

2.6390 

.78827 

.40938 

343 

Probabilistic 

1 

2.6389 

.77386 

.45011 

244 

Probabilistic 

2 

2.6379 

.78972 

.40427 

56 

every  case  the  optimization  results  were  essentially  the  same 
for  LGRAD  =  0  and  1.  The  composite  driveshaft  design  in 
Chapter  V  provides  an  excellent  example  showing  the  function 
evaluation  reductions  when  performing  a  relatively  complex 
design.  Table  IV  demonstrates  the  results  of  reducing  the 
frequency  of  constrained  variable  gradient  calculations 
when  designing  the  composite  driveshaft  for  minimum  weight, 
using  the  normal  distribution  model.  For  this  particular 
problem,  calculation  of  constrained  variable  gradients  only 
at  the  beginning  of  optimization  was  insufficient  to  obtain 
results  near  the  optimum. 


TABLE  IV 

REDUCTION  OF  FUNCTION  EVALUATIONS:  DRIVESHAFT  DESIGN 


LGRAD  WEIGHT  NUMBER  OF  FUNCTION 

LBS .  EVALUATIONS 


0 

1 

2 


9.0203 

9.0205 

27.216 


4011 

1344 
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VI.  NUMERICAL  EXAMPLES 


The  following  examples  are  presented  to  demonstrate  the 
capabilities  of  the  computer  subroutines  developed  in  this 
investigation.  In  both  cases  an  existing  ANALIZ  subroutine 
was  used  after  minor  modifications.  For  these  two  shaft 
designs  the  subroutine  from  [Ref.  6] ,  was  used.  The  ANALIZ 
subroutine  was  altered  to  include  the  required  probability 
variables  in  the  global  common  block  and  to  remove  the 
factor  of  safety  calculations  which  were  originally  used  as 
constraints . 

In  each  example  a  similar  deterministic  design  was  per¬ 
formed  to  give  the  reader  a  better  feel  for  the  results. 

The  deterministic  designs  are  not  intended  to  duplicate  the 
probabilistic  design  results. 

The  mean  values  used  in  the  designs  are  from  [Ref.  6] . 
The  coefficients  of  variation  are  reasonable  assumed  values 
used  for  demonstration  purposes. 

A.  STEEL  DRIVESHAFT 

Design,  for  minimum  weight,  a  steel  driveshaft  to  trans¬ 
mit  150.0  horsepower  at  300.0  RPM.  The  loading  is  presented 
in  Figure  6.1.  The  shaft  is  designed  against  failure  in 
strain,  in  torsional  buckling,  vibration  frequency  and  a 
maximum  deflection  of  0.05  inches. 

The  design  variables  are  the  shaft  thickness  and  the 
shaft  inside  diameter.  The  thickness  is  to  be  maintained 
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between  0.01  and  2.0  inches  and  the  inside  diameter  between 
0.5  inches  and  5.0  inches.  The  initial  design  is  an  inside 
diameter  of  2.0  inches  and  a  thickness  of  1.0  inch.  The 
remaining  data  for  the  analysis  is  presented  in  Table  V. 


TABLE  V 

INPUT  DATA  FOR  STEEL  DRIVESHAFT 


MEAN 

COEFFICIENT 

VALUE 

OF  VARIATION 

HORSEPOWER 

150. 

HP 

.08 

SPEED 

300  . 

RPM 

.08 

AXIAL  LOAD 

50. 

LB. 

.18 

LENGTH 

120  . 

IN. 

.015 

YOUNGS  MODULUS 

30,000,000 

PSI 

.038 

SHEAR  MODULUS 

11,538,000 

PSI 

.042 

POISSONS  RATIO 

.3 

.026 

SPECIFIC  WEIGHT 

.282 

LB/CU.  IN. 

.01 

STRAIN  LIMIT 

.001 

IN. /IN. 

.01 

MAX  DEFLECTION 

.05 

IN. 

.05 

THICKNESS 

.015 

INSIDE  DIAMETER 

.015 

FAILURE  DUE  TO: 

TORSIONAL  BUCKLING 

.21 

FREQUENCY 

.06 

Allowing  for  maximum  probabilities  of  failure  of  0.01  in 
strain,  torsional  buckling,  deflection  and  vibration  fre¬ 
quency;  a  minimum  weight  of  54.15  lbs  is  achieved  for  the 
normal  model  and  52.88  lbs  for  the  lognormal  model.  The 
initial  design  summary  is  presented  in  Table  VI  and  the 
final  designs  summaries  and  COPES  optimization  results  in 
Tables  VII  through  XII. 

Designing  for  a  factor  of  safety  against  failure  of  2.0 
for  strain,  torsional  buckling,  vibration  frequency,  and 


TABLE  VI 


STEBL  DRIVESHAFT  INITIAL  DESIGN  SUMHARY 


STS  EL  DRIVESHAFT  OUTPUT 


NUMBER  OF  PLYS  =  1 

NUMBER  OF  MATERIAL  TYPES  *  1 

NUMBER  OF  LOAD  CONDITIONS  =  1 

ECCENTRICITY  =  3.0 


DIMENSIONS 

PLY  THICKNESS  PERCENT  DIAMETER  THETA 

INSIDE  0.20000E+01 

1  0.  10000E+01  100.00  3.40003S+01  0.0 


STIFFNESS 
AE  =  0.31071 S+09 
El  =  0.38838E+09 
GJ  =  0.2718  7E+09 


LOADS  : 

L.C.  T  a  F  H.P.  3PM 

1  0.3151 3E+05  0.0  0.50000E+02  0.15000E+03  0.3000QE+03 


LOAD  CONDITION  1 

PLY  EPL  S.F.  SPT  S.F.  5PLT  S. F. 

1  0.16092E-06  100.00  0.0  100.00  0.23182E-03  5.61 

CRITICAL  SPEED  =  0.1554  3S+0’4 

MAXIMUM  DEFLECTION  =  3.19195E-G1 
WEIGHT  =  0 .31893E+03 

VOLUME  =0.1 13 10E+0U 
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TABLE  VII 


STEEL  DRIVESHAFT  COPES  STANDARD  OUTPUT:  NORMAL 


OPTIMIZATION  RESULTS 


OBJECTIVE  FUNCTION 
GLOBAL  LOCATION  3 


FUNCTION  VALUE  0.54154E+02 


DESIGN  VARIABLES 

D.  V.  GLOBAL 
ID  NO.  VAR.  NO. 
1  1  11 

2  2  2 


LONER 
30U  ND 

0.1030  OE-3  1 
0.5000  0E+00 


VALUE 

3. 99883E-01 
0. 50000E+01 


UPPER 

BOUND 

0 . 200  00E+0 1 
0.50000E+01 


DESIGN  CONSTRAINTS 

DETERMI NIS TIC 


ID 

1 

3 

5 

7 

8 
9 


GLOBAL 
VAR.  NO. 

51 

52 

53 

54 

56 

57 


LOWER 

BOUND 

■0.  10000  E-0 
■0.  10  000  E-0 
•0.  10000  E-0 
0. 31 513Ef 3 
0.30 000 £+0 


-0.  11000E+16 


MEAN 

VALUE 

3.9477UE-06 

0.0 

0.6821 3S-03 
0.  2434  0E  +  06 
0. 2506  9E«-04 
0.72076E-02 


U  P°E  R 
BOUND 

0.  10000E-02 
0. 100 00 E-0 2 
0. 10000E-02 
0.1  1 000  E*  1 6 
0.1  1030E+16 
G.50000E-0 1 
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TABLE  VIII 


STEEL  DRIVESHAFT  COPES  PROBABILISTIC  OUTPUT:  NORMAL 


ID 

1 

3 

5 

7 

8 
9 


PROBA3I LISTIC 


MEAN 

VALUE 

0.94774E-06 

0.0 

0.68213E-03 
0.24340E+06 
0. 2506  8E+  04 
0.72076E-02 


STANDARD 
DEVIATION 
0. 17535-06 
0.0 

0. 6535E-04 
0. 1364E  +  05 
0. 9634E+02 
0.57073-03 


COEFFICIENT 
OF  VARIATION 
0. 1753E-01 
0.0 

0. 9581E-01 
0. 5606E-01 
0. 38U3E-01 
0.791 8E-01 


PROBABILITY  OF  FAILURE 
AT  LONER  BOUND 


PROBABILITY  OF  FAILURE 
AT  UPPER  BOUND 


D 


ALLOWED 
0. 10 OOE-O 1 
0. 1C  OOE-O  1 
0. 10  OOE-O  1 
0. 10  OOE-O  1 
0. 10  OOE-O  1 
0. 10  OOE+O  1 


CALCULATED 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


0.  1  0  00E-( 
0.  1 0  00E-( 
0.  1  0  00E  +  ( 
0.  10005  +  ! 
0.  1 0  00E-( 


CALCULATED 
0.0 
0.  0 

0.  1 0  00E-01 
0.  0 
0.  0 
0.0 
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TABLE  II 


\ 

i 

I 

i 

I 


STEEL  DRIVESHAFT  DESIGN  SUMMARY:  NORMAL 


STEEL  DRIVESHAFT  OUTPUT 


NUMBER  OF  PLYS  =  1 

NUMBER  OF  MATERIAL  TIP ES  =  1 

NUMBER  OF  LOAD  CONDITIONS  =  1 

ECCENTRICITY  =  0.0 


DIMENSIONS 

PLY  THICKNESS  PERCENT  DIAMETER  THETA 

INSIDE  0.50000E  +  0  I 

1  0. 99883E-0 1  100.00  0.51993E+01  0.0 


STIFFNESS 
AE  =  0.52757E+08 
El  =  0.171 58E+09 
GJ  =  0.12011E+09 


LOADS  : 

L.C.  T  a 

1  0. 31 5 13E+Q5  0.0 


0.50  0OOE+O2 


H.P. 

0. 1 5000E+03 


LOAD  CONDITION  1 

PLY  EPL  S.F.  EPT  S.F. 

1  0.94774E-06  100.00  0.0  100.00 


CRITICAL  SPEED 
MAXIMUM  DEFLECTION 
HEIGHT 
VOLUME 


=  0 .5U154E+02 
=  0.192  04E+03 


0.68213E-0 


) 


i 


* 

i 
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RPM 

. 300  00E+0  3 


S.F. 
3  1.9  1 


TABLE  X 
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T ABLE  XI 


STEEL  DRIVESHAPT  COPES  PROBABILISTIC  OUTPUT:  LOGNORMAL 


PROBABILISTIC 


ID 

1 

3 

5 

7 

8 
9 


HEAR 

VALUE 

0.  97054E-06 

0.0 

0. 69856E-03 
0. 2296 1E+06 
0.  25056E  +  04 
0.  72144E-02 


STANDARD 
DEVIATION 
0. 1796B-06 
0.0 

0.6692E-04 
0. 1283E+05 
0.9630E+02 
0. 57 1 3 E- 03 


COEFFICIENT 
OF  VARIATION 
0.  1796E-0 1 
0.0 

0. 9580E-01 
0.561  IE-01 
0 . 3843E-0 1 
0.791 8E-01 


PROBABILITY  OF  FAILURE 
AT  LOWER  BOUND 


ID  ALLOWED 

0.  1000E-01 
0.  10  OOE-0 1 
0. 1000E-01 
0.  1000E-01 
0.  1000E-01 
0.  1000E+01 


CALCULATED 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


PROBABILITY  OF  FAILURE 
AT  UPPER  BOUND 


ALLOWED 
0. 1000E-01 
0.  1000E-01 
0.  1 0  00E-01 
0. 1 0  OQE  +  0 1 
0. 1000S+G1 
0.  1000E-01 


CALCULATED 

0.0 

0.0 

0. 9993E-02 
0.  0 
0.0 
0.0 
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TABLE  XII 


STEEL  DRIVES HAPT  DESIGN  SUMMARY:  LOGNORBAL 


STEEL  DRIVESHAPT  OUTPUT 


NUMBER  OF  PLYS  =  1 

NUMBER  OP  MATERIAL  TYPES  =  1 

NUMBER  OF  LOAD  CONDITIONS  =  1 

ECCENTRICITY  =  3.0 


DIMENSIONS 

PLY  THICKNESS  PERCENT  DIAMETER  THETA 
INSIDE  0.50  000E  +  01 

1  0.  975 8 1 E-0 1  103.00  0.51952E+01  0.0 


STIFFNESS 
AE  =  0. 515182+08 
El  =  0.  1 674  0E+  09 
GJ  =  0.11718E+09 


LOA  DS  : 

L  " 

i~*0. 31513  E+ 05  0.0 


0.53  000E  +  02 


H.  P. 

0. 15000E+03 


RPM 

0.30000E+03 


LOAD  CONDITION  1 
PLY  EPL  S.F. 

1  0. 97054 E-06  100.00 

CRITICAL  SPEED 
MAXIMUM  DEFLECTION 
HEIGHT 
VOLUME 


0.0 


PT 


S.  F. 

100. 00 


=  0 . 25056E+04 
=  0 .721  ‘ 


=  0 
=  3 


44E-02 
52882E+02 
18752E+03 


EPLT  S.F. 

0. 69856E-03  1.86 
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deflection  an  optimum  of  73.78  lbs  is  obtained.  The  COPES 
optimization  results  and  final  design  summary  are  contained 
in  Tables  XIII  and  XIV. 

B.  COMPOSITE  DRIVESHAFT 

Design,  for  minimum  weight  a  120  in.  long,  four  ply 
graphite  epoxy  driveshaft  to  transmit  150.0  horsepower  at 
300  RPM.  Figure  6.1  shows  the  shaft  loading  condition. 

The  shaft  is  designed  against  failure  in  transverse,  logitu- 
dinal,  and  shear  strain;  torsional  buckling,  vibration  fre¬ 
quency  and  a  maximum  deflection  of  0.05  inches.  Table  XV 
contains  the  initial  design. 

The  design  variables  are  the  thickness  of  each  ply,  the 
orientation  of  each  ply  and  the  inside  diameter.  Addition¬ 
ally,  the  thickness  of  ply  2  and  ply  3  must  remain  equal  to 
each  other  and  the  orientation  of  ply  3  equals  the  negative 
of  ply  2.  The  following  additional  constraints  are  placed 
on  the  design  variables;  The  inside  diameter  must  remain 
between  0.5  inches  and  5.0  inches  and  the  thickness  of  any 
given  ply  between  0.01  and  0.5  inches.  Ply  2  is  allowed  to 
vary  between  0.0  and  90.0  degrees.  The  orientation  of  plys 
1  and  4  remain  constant.  The  remaining  input  data  is  pre¬ 
sented  in  Table  XVI. 

Allowing  for  maximum  probabilities  of  failure  of  0.01  in 
strain,  torsional  buckling,  deflection,  and  vibration  fre¬ 
quency;  an  optimum  weight  of  9.02  lbs  is  achieved  using  the 
normal  distribution  model  and  10.28  lbs  for  the  lognormal 
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TABLE  XIII 


STEEL  DHIVESHAFT  COPES  OUTPUT :  DETERMINISTIC 


OPTIMIZATION  RESULTS 


OBJECTIVE  FUNCTION 

GLOBAL  LOCATION  3  FUNCTION  VALUE  0.73783E+02 


DESIGN  VARIABLES 

D.  V.  GLOBAL 
ID  NO.  VAR.  NO. 
1  1  11 


LONER 

30UND 

0.1030  OE-O 1 
0.5030  OE+Q  0 


VALUE 

3.  13515E+00 
0. 5000QE+01 


UPPER 

BOUND 

0. 20000E+01 
0.50000E+01 


DESIGN  CONSTRAINTS 

DETERMI  NISTIC 


GLOBAL 

LONER 

MEAN 

UPPER 

ID 

VAR.  NO. 

BOUN  D 

VALUE 

BOUND 

1 

51 

-0.50000E-03 

0.5956  IE-06 

Q.50000E-03 

3 

52 

-0.  50000 E-03 

0.3 

0 . 5000  OE-O  3 

5 

53 

-0. 53000  E-03 

0. 5003  5S-03 

0 . 50000E-03 

7 

54 

0  .63025  E+05 

0. 518382+06 

3. 11000E+16 

8 

56 

0.6J303  E+03 

0.25247E+04 

0. 11000E+16 

9 

57 

-0 . 1 1 000  E+  16 

0.7 1 044E-02 

0 . 20 00 0E*0 1 
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TABLE  XIV 


STEEL  DHIVESHAPT  DESIGN  SUMMARY:  DETERMINISTIC 


STEEL  DRIVESHAFT  OUTPUT 


NUMBER  OF  PLYS 

NUMBER  OF  MATERIAL  TYPES  = 

NUMBER  OF  LOAD  CONDITIONS  = 

ECCENTRICITY 


DIMENSIONS 

PLY  THICKNESS  PERCENT  DIAMETER 
INSIDE  0.50000E  +  01 
1  0. 1 35 15E+00  100.00  Q.52703E+01 


THETA 


STIFFNESS 
AE  =  0.7 1879E+08 
El  =  0.237  09E+09 
GJ  =  0.165  97E+09 


LOADS: 

Lr'0.315l3E+05 


0.50  000E  +  02 


10 

0. 1 5000E+03 


RPM 

.300C0E+03 


LOAD  CONDITION  1 

PLY  EPL  S.F.  EPT  S.F. 

1  0. 69561 E-06  100.00  0.0  100.00 

CRITICAL  SPEED  =  0.25247E  +  04 

MAXIMUM  DEFLECTION  =  0.71044E-02 

HEIGHT  ■  0.73783E+02 

VOLUME  =  0.26164E+03 


EPLT  S.F. 

0 . 50035E-03  2.60 
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TABLE  XV 


COMPOSITE  DRIVESHAFT  INITIAL  DESIGN 


THICKNESS 

IN. 


ORIENTATION 

DEGREES 


PLY  1 
PLY  2 
PLY  3 
PLY  4 


25 

0.0 

25 

20.0 

25 

-20.0 

25 

90.0 

INSIDE  DIAMETER  2.0  IN. 


TABLE  XVI 

INPUT  DATA  FOR  COMPOSITE  DRIVESHAFT 


MEAN  COEFFICIENT 

VALUE  OF  VARIATION 


HORSEPOWER 

SPEED 

LENGTH 

AXIAL  LOAD 

LONGITUDINAL  MODULUS 

TRANSVERSE  MODULUS 

SHEAR  MODULUS 

MAJOR  POSSIONS  RATIO 

STRAIN  LIMITS: 

LONGITUDINAL  COMPRESSIVE 
LONGITUDINAL  TENSILE 
TRANSVERSE  COMPRESSIVE 
TRANSVERSE  TENSILE 
MAX  SHEAR 
SPECIFIC  WEIGHT 
MAX  DEFLECTION 
PLY  ORIENTATION 
PLY  THICKNESS 
INSIDE  DIAMETER 
FAILURE  DUE  TO: 

TORSIONAL  BUCKLING 
FREQUENCY 


150. 

HP 

.08 

300  . 

RPM 

.08 

120. 

IN. 

.015 

50. 

LB. 

.18 

21,000,000 

PS  I 

.14 

1,700,000 

PS  I 

.14 

650,000 

PS  I 

.14 

.21 

.08 

-.00857 

IN. /IN. 

.15 

.00857 

IN. /IN. 

.15 

-.0176 

IN. /IN. 

.15 

.00471 

IN. /IN. 

.15 

.0184 

IN. /IN. 

.15 

.056 

LB/CU.  IN. 

.05 

.05 

IN. 

.05 

.05 

.015 

.015 

.20 

.06 
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model.  The  ANALIZ  subroutine  design  summary  for  the  initial 
design  is  presented  in  Table  XVII.  Tables  XVIII  through 

XXIII  contain  the  design  summaries  and  final  COPES  optimiza¬ 
tion  results. 

Designing  for  factors  of  safety  against  failure  of  2.0 
for  strain,  torsional  buckling,  deflection,  and  vibration 
frequency  an  optimum  weight  of  10.23  lbs  is  obtained.  Tables 

XXIV  and  XXV  present  the  deterministic  design  results. 
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TABLE  XVII 


COMPOSITE  OBI TESHAET  INITIAL  DESIGN  SUMMARY 


COMPOSITE  DRI VESHAFT  OUTPUT 


NUMBER  OF  PLYS  =  4 

NUMBER  OF  MATERIAL  TYPES  =  1 

NUMBER  OF  LOAD  CONDITIONS  =  1 

ECCENTRICITY  =  3.0 


DIMENSIONS 

PLY  THICKNESS  PERCENT  DIAMETER 
INSIDE  0.20  000E  +  0 1 

1  0. 25000E+00  25.00  0.25000E+01 

2  0. 25000E+00  25.00  0.30000E+01 

3  0. 250 00E+0 0  25.00  0.35000E+01 

4  0.  250 00E+00  25.00  0.4000DE+01 


THETA 

0.0 

0. 20000E+02 
-0. 200QQE+02 
0 . 90000E+02 


STIFFNESS 
AE  =  0.1 2143E+  09 
El  *  0. 12426E+09 
GJ  =  0 . 3722  5E*  08 


OADS: 


L.C 

T 

M 

F 

H. 

P.  RPM 

1 

0.31513E+05 

0.  0 

0.50000E  +  02 

0 .  1  5000S+03  0 . 300  00E  +  0  3 

LOAD  CONDITION 

1 

PLY 

EPL 

S.  F. 

EPT 

S.  F. 

EPLT 

S.P. 

1 

0. 41 176E-06 

100. 00 

0.0 

100. 00 

4  3.13 

0. 10582E-02 

17.39 

2 

C. 40847E-03 

20.98 

-0.40  806B-03 

0.97246E -03 

18.92 

3  - 

0. 47576 S-03 

18.  01 

0.47617E-03 

9.89 

0. 11351E-02 

16.21 

4 

0. 53 159E-09 

100. 00 

0.41  123E-06 

100.00 

-0.  16531  E- 02 

10.87 

CRITICAL  SPEED 

=  0 . 19725E+ 

04 

MAXIMUM  DEFLECTION  =  0.11746E- 

01 

WEIGHT 

VOLUME 

=  0 . 63334E+02 
=  0.1 1310E+04 
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T4BL2  IVIII 


COMPOSITE  DRI7ESH1FT  COPES  STANDARD  OUTPUT:  NORMAL 


OPTIMIZATION  RESOLTS 


OBJECTIVE  FUNCTION 

GLOBAL  LOCATION  3  FUNCTION  VALUE  0.90203E+01 


DESIGN  VARIABLES 


D.  V. 

GLOBAL 

LO  HER 

UPPER 

ID 

NO. 

VAR.  NO. 

BOUND 

VALUE 

BOUND 

1 

1 

11 

0. 1030  OE-O 1 

0.  10000E-01 

0.50  0  00E+00 

2 

2 

12 

0.  1000 OE-O  1 

3.  1C212E+00 

0.50000E+00 

3 

2 

13 

0.  10000E-01 

0.  10212E+00 

0.500 OOE+QO 

a 

3 

14 

0. 1000  OE-O  1 

0.  10000E-01 

0 .5  00  00E  +  00 

5 

a 

32 

0.0 

0. 34968E+02 

0 . 900  90E  +  02 

6 

4 

33 

0.0 

-0. 3U368E+02 

-0.90000E+Q2 

7 

5 

2 

0.50000E+00 

0.  16311E+01 

0.50000E+01 

DESIGN 

CONSTRAINTS 

DETERSI  NISTIC 

GLOBAL 

LOHER 

MEAN 

UPPER 

ID 

VAR.  NO. 

BOUND 

VALUE 

BOUND 

1 

51 

-0. 85700  E-02 

0.  35785E-05 

0.85700E-02 

3 

52 

-0. 17600  E-01 

0.0 

0.47 100E-02 

5 

53 

-0 .18400  E-01 

0.47679E-02 

0.  18400E-01 

7 

54 

-0  .85  700  E-02 

0.251 05E-02 

0 . 85700E-02 

9 

55 

-0.17600  E-01 

-0.25069E-02 

0.471 00E-02 

11 

56 

-0.19400 E-01 

0.  13288E-02 

0.  18400E-01 

13 

57 

-0.85700  E-02 

-0. 27745E-02 

0.857  OOE-02 

15 

58 

-0 .17600E-01 

0.2  778  IE-02 

0.47  100E-02 

17 

59 

-0.18400 E-01 

0.203 19E-02 

0.  1840QE-01 

19 

60 

-0 .85700E-02 

0.  18741E-08 

0 . 857  3QE-02 

21 

6  1 

-0 . 17600  E-01 

0.  35766E-05 

0.471 Q0E-02 

23 

62 

-0 . 18400  E-01 

-0.  59688E-02 

0.  18400 E-01 

25 

63 

0 . 31  51  3  E+05 

0.55 160E+05 

0. 11030E+16 

26 

65 

0 .30000E+03 

0.  1 17 15E+04 

0.  11000E+16 

27 

66 

-0 .11000  E+16 

0. 34815E-01 

0.50000E-01 
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r&BLB  ZII 


COMPOSITE  DRI VESHAFT  COPES  PROBABILISTIC  OOTPOT:  NORMAL 


PSOBABI  LI  STIC 


MEAN 

STANDARD 

COEFFICIENT 

ID 

V  AIDE 

DEVIATION 

OF  VARIATION 

1 

0 . 35785E- 05 

0.4929E-36 

0. 4929E-01 

3 

0.0 

0.0 

0.0 

5 

0.47679E-02 

0.7317E-03 

0. 1535E+00 

7 

0.25105E-02 

0. 38385-03 

0.  1529E+00 

9 

-0.25069E-02 

0. 38345-03 

3.  1530E+C0 

1 1 

0. 18288 E- 02 

0. 4363E-03 

0. 2386E+00 

13 

-0 . 27745E- 02 

0.422 15-03 

0.  1521 2*00 

15 

0 . 2778 1 E- 02 

0.4225E-03 

0.  1521E+00 

17 

0.2031 9E- 02 

0.4938E-03 

0. 2430E+00 

19 

0. 18741E-08 

0.46875-03 

0. 1000E+02 

21 

0. 35766E- 05 

0.4687E-93 

0.  1000E+02 

23 

-0 . 59688  E- 02 

0. 9181E-03 

0.  1538E  +  00 

25 

0.5516  0E+  05 

0. 78492*04 

3  .  1423E+00 

26 

0. 1 1715E+04 

0. 9463E+02 

0. 8078E-01 

27 

0.34815E-01 

0.6032E-02 

3.  1733E+00 

PROBABILITY 

OF  FAILURE 

PROBABILITY 

OF  FAILURE 

AT  LOWER 

BOUND 

AT  UPPER 

BOUND 

CALgULAT ED 

C^LgULATED 

3 

0. 1000E-01 

o.'o 

0.  10005-01 

ol  0 

5 

0.  1000E-01 

0.0 

0.  1 3  OOE-0 1 

0. 9537E-06 

7 

0. 1000E-01 

0.0 

0. 1 0  OOE-0 1 

0. 3338E-05 

9 

0. 1000E-01 

0  .0 

0.  1 0  OOE-0 1 

0.0 

11 

0. 1000E-01 

0.0 

0. 1 0  OOE-O 1 

0.  0 

13 

0. 1000E-01 

0.95  37E-05 

0.  1  0  OOE-O 1 

0.  0 

15 

0. 10005-01 

0.0 

0.  1 0  OOE-0 1 

0. 9468E-02 

17 

0.  1000E-01 

0  .0 

0.  1000E-01 

3.0 

19 

0. 10  OOE-O  1 

0.0 

0. 1 0  OOE-0 1 

0.  0 

21 

0.  1000E-0  1 

0.0 

0. 1 0  OOE-0 1 

0.  0 

23 

0. 1000E-0 1 

0 . 10  01 E -04 

0.  1 0  OOE-0 1 

0.  0 

25 

0. 10  OOE-O  1 

0.94105-02 

0.  1300E  +  01 

0.  0 

26 

0. 10  OOE-O  1 

0.0 

0. 1 0  OOE  +  0 1 

0.  0 

27 

0.  10  00E  +  0  1 

0.0 

0.  1 0  OOE-O 1 

0.  1902E-01 
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TABLE  XX 


COMPOSITE  OSI VESHAFT  DESIGN  SUMMARY:  NORMAL 


COMPOSITE  DRI VESHAFT  OUTPUT 


NUMBER  OP  PLYS  =  4 

NUMBER  OF  MATERIAL  TYPES  =  1 

NUMBER  OF  LOAD  CONDITIONS  =  1 

ECCENTRICITY  =  3.0 


DIMENSIONS 


PLY 

1 

2 

3 

4 


0. 

0. 

0. 

0. 


THICKNESS  PERCENT  DIAMETER 
INSIDE  0. 15  81 1E+01 
10000E-01  4.45  0.17011E+01 

102 12E+00  45.54  0.19054E+01 

102 12E+00  45.54  0.21095E+01 


100  OOE-0 1 


STIFFNESS 
AS  =  0.  139  72E+  08 
El  =  0.631 16E+07 
GJ  =  0 . 562 1 7E+07 


4 .45 


3 . 2 1 295  E*  0  1 


THETA 

0.0 

0. 34  96  SE  +  02 
-0.34968E+02 
0. 93000E  +  02 


LOADS  : 

L  ^  T  M 

1*0.3151 3E+05  0.0 


LOAD  CONDITION  1 
PLY  SPL  S.F 

1  0.  35785E-05  100.  00 

2  0.25105E-0  2  3.  41 

3  -0.  277  45  E-Q  2  3. 09 

4  0.  18741 E-08  100.  00 


0.50  0000E  +  02  3 


H  ® 

1  5000E+03 


RPM 

30000E+33 


SPT  S.F. 

0.0  130.00 

-0.25  069E-02  7.  02 

0.27  781E-02  1.70 

0.35  766E-05  100.  00 


EPLT  S.F. 

0.47579E-02  3.86 

0.  18288E-02  10.06 
0.203  1 9E-0  2  9.06 

-0.59688E-02  3.08 


CRITICAL  SPEED 
MAXIMUM  DEFLECTION 
HEIGHT 
VOLUME 


=  3. 11715  E*  04 
=  3 . 3  48  15  E- 01 
=  0 .902  03E+01 
=  0  .  16*  i8E*03 
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TABLE  XXI 

COMPOSITE  D8I7ESHAPT  COPES  STANDARD  OUTPUT:  LOGNORMAL 


OPTIMIZATION  RESULTS 
OBJECTIVE  FUNCTION 

GLOBAL  LOCATION  3  PUNCTION  VALUE  0.10281E+Q2 


DESIGN  VARIABLES 


D.  V. 

GLOBAL 

LOW  ER 

UPPER 

ID 

NO. 

VAR.  NO. 

30U  ND 

VALUE 

30UND 

1 

1 

11 

0.  ICOOOE-O  1 

0.  19474E-01 

0 . 5  00  OOE  +  OO 

2 

2 

12 

0. 1000  OE-O 1 

0.  92916E-01 

0 .500  OQE  +  QO 

3 

2 

13 

0.  130002-0  1 

0. 92916E-01 

0 .5  00  00E+  00 

4 

3 

14 

0.1  00  0  OS-O  1 

0.  10000E-0 1 

0.50000S+00 

5 

4 

32 

0.0 

0.  41392E+02 

0.90000E+02 

6 

4 

33 

0.0 

-0.  41  392E+02 

-0.90000E+02 

7 

5 

2 

0.5000  0E+00 

0.  20465E  +  0  1 

0 . 5  00O0E+01 

DESIGN  CONSTRAINTS 

DETERMI NISTIC 


GLOBAL 

LOWER 

MEAN 

UPPER 

ID 

VAR.  NO. 

BOON  D 

VALUE 

BOUND 

1 

51 

-0.85700  E-02 

0.  37986E-05 

0.85700E-02 

3 

52 

-G. 17600 E- 01 

0.0 

0 . 47 100E-02 

5 

53 

-0 .18400  E-0 1 

0.34484E-02 

0.  18400E-01 

7 

54 

-0.85  700  E-02 

0.  18651 E-02 

0 . 85  7 OOE-O  2 

9 

55 

-0 . 17600  E-01 

-0.  1 3613E-02 

0.4710  0E- 02 

1  1 

5o 

-0 . 18U00 E-01 

0. 46804S-03 

0.  18400E-01 

13 

57 

-0.85  700  E-02 

-0. 20  1  32S-02 

0.85700E-02 

15 

58 

-0 . 17600  E-01 

0.20170E-02 

0.47  100 E-0 2 

17 

59 

-0 . 18400  E-01 

3.51U18S-03 

0.  18400  E-01 

19 

80 

-0 .85700  £-02 

0.  1  286  1E-08 

0. 85700E-02 

21 

6  1 

-0 . 17600  E-01 

0.37974E-05 

0. 47 100E-02 

23 

62 

-0 . 18U00  E-01 

-0 . 4 0960 E-02 

0.  18400E-01 

25 

63 

0 .31  513  E  +  05 

0. 549  75E  +  05 

0.11  OQOE+16 

26 

65 

0 .30000  E+03 

0.  12564S  +  0U 

0.  1  1000E+ 16 

27 

5  6 

-0 . 11000  2+  1  6 

0. 29994E-01 

0. 50000E-01 

TABLE  XXII 


CO MPOSITE  DRIVESHAFT  COPES  PROBABILISTIC  OUTPUT:  LOGNORMAL 


PROBABILISTIC 


MEAN 

STANDARD 

COEFFICIENT 

ID 

VALUE 

DEVIATION 

OF  VARIATION 

1 

0. 37986E-05 

0.5215E-06 

0. 52 15E-01 

3 

0.  0 

0.0 

0.0 

5 

0. 34484E -02 

0.5257E-03 

0. 1524E+00 

7 

0.  1865  IE-02 

0.  2842E-03 

0. 1 524E+00 

9 

-0.  1861  3E-02 

0. 2838E-03 

0.  15  25E+00 

11 

0.  46804E-03 

0 . 28  1 5  E-  0  3 

0.601 5E+00 

13 

-0.  20  1  3  2E-02 

0.3069  E-03 

0.  1524E+00 

15 

0. 2017QE-02 

0.3073  E- 03 

0. 1 524E+00 

17 

0. 5141 8E-03 

0. 3062E-03 

0 . 5956E+00 

19 

0-  1236  IE-08 

0.3215 E-03 

0. 1000E+02 

21 

0- 37974E -05 

0.3216  E-03 

0.  1000E+02 

23 

-0. 40960E-02 

0. 6256  E- 03 

0.  1527E+00 

25 

0.  54975E  +  05 

0.7813E+04 

0.  1421E+00 

26 

0.  12564E+04 

0. 1012E+03 

0. 8058E-01 

27 

0. 29994E-01 

0.513  5  E-0  2 

0. 1712E+00 

PROBABILITY 

OF  FAILURE 

PROBAB 

TT 
M.  Xj 

ITY  OF  FAILURE 

AT  LOWER 

BOUND 

AT 

UPPER  BOUND 

ID 

ALLOWED 

CALCULATED 

ALLOWED 

CALCULATED 

1 

0.  1000E-01 

0.0 

0.  1000E- 

01 

0.0 

3 

0.  1000E-01 

0.0 

0.  1000E- 

01 

0.0 

5 

0.  1000E-01 

0  .0 

0.  1000E- 

01 

0.0 

7 

0.  1000E-01 

0  .0 

0. 1000E- 

01 

0.0 

9 

0. 10  COE-O  1 

0.0 

0. 1000E- 

01 

0.0 

1  1 

0. 1000E-0 1 

0.0 

0.  1000E- 

01 

0.0 

13 

0.  1000c  '01 

0.0 

0. 1000E- 

01 

0.0 

15 

0. 10  OOE-O  1 

c.o 

0. 1Q00E- 

01 

0.3333E-04 

17 

0.  10  OOE-O  1 

C  .0 

0. 1000E- 

01 

0.0 

19 

0. 1000E-01 

0.0 

0.  1000E- 

01 

0.531  3E-02 

21 

0. 10  OOE-O 1 

0.0 

0.  1000E- 

01 

0.9957E-02 

23 

0. 10  OOE-O  1 

0.0 

0.  1000E- 

01 

0.0 

25 

0. 10  OOE-O  1 

0 . 1000E -01 

0. 1000E+01 

0.0 

26 

0.  1000E-0  1 

0.0 

0.  1000E  +  01 

0.0 

27 

0.  10  00E+-0  1 

0  .0 

0.  1000E- 

01 

0. 1 54  4  E-0  2 

68 


TABLE  XXIII 


COS  POSITS  DRIVESHAFT  DESIGN  SUMMARY:  LOGNORMAL 


COMPOSITE  DRIVESHAFT  OUTPUT 


NUMBER  OF  PLYS  =  4 

NUMBER  OF  MATERIAL  TYPES  =  1 

NUMBER  OF  LOAD  CONDITIONS  =  1 

ECCENTRICITY  =  0.0 


DIMENSIONS 

PLY  THICKNESS  PERCENT  DIAMETER 
INSIDE  0.20465E+01 
0.  19474E-01  9.04  0.20855E+Q1 

0.  92916E-01  43.16  0.22713E+01 

0.  92916E-01  43.16  0 . 24571 E+01 


1 

2 

3 

4 


0.  100  OOE-Q  1 


4.54 


0 .24771 E+  01 


THETA 

0.0 

0.41 39  2E+0  2 
-0.41 39  2E+02 
0. 90000E+02 


STIFFNESS 
AE  =  0.131 63E+08 
El  =  0.82517E+07 
GJ  =  0.952  89E+  07 


L 


LOADS  : 
0.31513E+05 


M  F  H.P.  RPM 

0.0  0 . 5 0 00 00 E  +  02  2.  15000S+03  0.30000E+03 


LOAD  CONDITION  1 

PLY  E  PL  S.F.  SPT  S.F. 

1  0.  37986E-05  100.  00  0.0  100.00 

2  0. 18651 E-02  4.59  -0.18613E-02  9.46 

3  -0.  20  1  32E-02  4.26  0.20  170E-02  2.34 

4  0.  12861  E-0  8  100.00  0.37974E-05  10  0.  00 

CRITICAL  SPEED  =  0.125643*04 

MAXIMUM  DEFLECTION  =  0.29994E-01 

WEIGHT  =  0  .  1028  1E+02 

VOLUME  =  0  .  1 83  5  9E  +  03 


EPLT 

0. 34484E-02 
0. 46804E-03 
0. 51418 E-03 
-0. 40960E- 02 


5. 34 
39.31 
35.79 
4.49 
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TABLE  XXI? 


COMPOSITE  DRIVESHAFT  COPES  OUTPUT:  DETERMINISTIC 


OPT IM IZATION  RESULTS 


OBJECTIVE  FUNCTION 

GLOBAL  LOCATION  3  FUNCTION  VALUE  0.  1023 1E  +  02 


DESIGN  VARIABLES 


ID 

D.  V. 
NO. 

GLOBAL 
VAR.  NO. 

LOW  ER 
BOUND 

7ALUE 

UPPER 

BOUND 

1 

1 

11 

0.10000E-0  1 

0.  1  OOOOE-0 1 

0.50000E+00 

2 

2 

12 

0.10000E-0  1 

0.  1 06  97  E+00 

0 . 50  3  00E+00 

3 

2 

13 

0.  10300E-01 

0.  10697E+00 

0.50000E+00 

4 

3 

14 

0.1030  OE-O  1 

3. 10000E-01 

0.50000E+00 

5 

4 

32 

0.0 

0.  30524E+02 

0.90000E+02 

6 

4 

33 

0.0 

-0.  30524 E+ 02 

-0. 90000E+02 

7 

5 

2 

0. 5030  OE+Q  0 

0.  13376E+01 

0.50000E+01 

DESIGN 

CONSTRAINTS 

DETERS I  NISTIC 

GLOBAL 

LOWES 

MEAN 

UPPEP 

ID 

VAR.  NO. 

BOUND 

VALUE 

BOUND 

1 

51 

-0.42850  E-02 

0.26961 E-05 

0. 42850E-02 

3 

52 

-0.88000S-02 

0.3 

3. 23550E-02 

5 

53 

-0.92000  E-02 

0.43686E-02 

0. 92000E-02 

7 

54 

-0  .42  850  E-02 

0.2 1334E-02 

3.42850E-02 

9 

55 

-0.88000 E-02 

-0.21307E-02 

0. 23550E-02 

11 

56 

-0.92 000 E-02 

0.23560E-02 

0. 92000E-02 

13 

57 

-0.42  850  E-02 

-0.234  96E-02 

0.  4285 0E-02 

15 

58 

-0 .88000  E-02 

0.  23522E-02 

0.23  5  50E-02 

17 

59 

-0.92000  E-02 

0.26043E-02 

0. 92000F-02 

19 

60 

-0 .42  850  E-02 

0.  17024E-08 

0 . 42850E-02 

21 

6  1 

-0.88000  E-02 

0. 25944E-05 

0. 23550E-02 

23 

62 

-0 . 92000  E-02 

-0. 54219E-02 

0. 92000E-02 

25 

63 

0 .63025  E+05 

0.54108E+05 

0. 11 000E+16 

26 

65 

0  .60000  E+G3 

0.  1  3824S  +  04 

0. 1 1009E+ 16 

27 

66 

-0 .11000  E*16 

3.24515E-01 

0.25000E-01 
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TABLE  XZ7 


COMPOSITE  DRIVES HAPT  DESIGN  SUMMARY:  DETERMINISTIC 


COMPOSITE  DRI VESHAFT  OUTPUT 


NUMBER  OF  PL  IS  =  4 

BOMBER  OF  MATERIAL  TYPES  =  1 

NUMBER  OF  LOAD  CONDITIONS  =  1 

ECCENTRICITY  =  3.0 


DIMENSIONS 

PLY  THICKNESS  PERCENT  DIAMETER 
INSIDE  0 . 1 8  376E  +  0 1 
0.  10000E-01  4.27  0.18576E+01 

0. 10697E+00  45.73  0.20715E+01 

0.  10697E+00  45.73  0. 22855 E+ 01 

4.27  0. 23055E+  01 


1 

2 

3 

4 


0.  10000E-0  1 


THETA 

0.0 

0. 30524E+02 
-0.30524E+02 
0. 90000E+02 


STIFFNESS 
AE  =  0.18546E+08 
El  =  0.99273E+07 
GJ  =  0.66999E+07 


LOADS  : 

L.C.  T  a 

1  0.3151 3E+05  0.0 


0.50  0000E+02  0 


H.P. 

15000E+03 


RPM 

0. 30000E+03 


LOAD  CONDITION 


PLY 

EPL 

S.F. 

EPT 

S.F. 

SPLT 

S.F. 

1 

0. 26961 E-05 

100. 00 

0.0 

100.00 

0.43586E-02 

4.  21 

2 

0. 21334  E-02 

4. 02 

-0.21 307E-02 

8.  26 

0. 23560E-02 

7.81 

3  - 

0.  23496  E-0  2 

3.65 

0.23522S-02 

2.00 

0.  26343E-02 

7.07 

4 

0.  17024  E-08 

100.00 

0.26  944E-05 

100. 00 

-0.54219E-02 

3.39 

CRITICAL  SPEED  =  0.13824E+04 
MAXIMUM  DEFLECTION  =  0.24515E-01 
WEIGHT  =  0 . 1023  1E+02 
VOLUME  =  0.13270S+03 
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VII.  CONCLUSIONS 


Numerical  optimization  using  probabilistic  design  tech¬ 
niques  provides  an  effective  method  for  designing  structures 
based  on  allowed  component  reliability. 

Existing  ANALIZ  subroutines  can  be  easily  modified  to 
perform  probabilistic  optimization  and  as  demonstrated, 
complex  designs  can  be  accomplished. 

The  major  drawback  to  the  program  developed  is  the  large 
increase  in  function  evaluations  as  compared  to  the  standard 
deterministic  optimization.  This  disadvantage  can  be 
partially  offset  by  reducing  the  frequency  in  which  the 
constrained  variable  gradients  are  calculated. 


i 

i 
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VIII.  RECOMMENDATIONS 


A  method  of  scaling  the  variables  should  be  investigated 
to  reduce  ill-conditioning  in  the  optimization  process. 

This  is  a  possible  area  to  pursue  in  order  to  alleviate 
problems  caused  by  very  small  coefficients  of  variation. 

Further  efforts  should  be  undertaken  to  accomplish  the 
reduction  of  function  evaluations  while  maintaining  suitable 
optimization  results.  The  altering  of  the  frequency  in  which 
constrained  variable  gradients  are  calculated  appears  to  be 
the  most  promising  approach. 

In  order  to  reduce  function  evaluations  it  is  recommended 
that  COPES  be  modified  to  allow  the  user  to  supply  mathemati¬ 
cal  expressions  for  the  partial  derivatives  of  the  constrained 
variables  via  the  ANALIZ  subroutine.  These  could  be  used 
instead  of  the  finite  different  gradients  calculated  by 
COPE19/COPE20 .  Modifi cat ions  to  allow  the  user  to  supply 
precalculated  gradients  would  also  be  helpful  to  this  end. 

The  sophisticated  user  can  accomplish  this  now  by  using  CONMIN 
directly  in  conjunction  with  C0PE19 ,  C0PE21,  and  COPE22. 

It  is  further  recommended  that  additional  work  be  under¬ 
taken  to  modify  COPES  to  allow  for  the  use  of  the  following 
options  with  probabilistic  design:  A)  two  variable  function 
space;  B)  sensitivity  analysis;  C)  optimum  sensitivity. 

An  additional  worthwhile  modification  would  be  to  provide 
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for  the  calculation  in  C0PE19  and  subsequent  output  of  the 
standard  deviation  of  the  objective  function  mean  value. 

This  investigation  dealt  only  with  component  reliability 
in  the  design  process .  Work  should  be  undertaken  to  ascer¬ 
tain  the  potential  for  the  inclusion  of  system  reliability 
considerations  in  the  optimization  process. 
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APPENDIX  A 


COPES  MANUAL  ADDENDUM 

This  appendix  is  intended  as  an  addendum  to  the  COPES 
Manual,  [Ref.  1]. 


I.  INTRODUCTION 

The  purpose  of  this  document  is  to  provide  user  instruc¬ 
tions  for  performing  numerical  optimization  using  probabilis¬ 
tic  design  techniques.  Four  subroutines  have  been  included 
in  COPES  (Control  Program  for  Engineering  Synthesis)  to 
provide  this  additional  design  tool. 

This  discussion  describes  the  capabilities  of  probabilis¬ 
tic  optimization  using  the  COPES/CONMIN  program.  A  simple 
design  example  is  first  presented  to  demonstrate  the  program 
capabilities.  Guidelines  are  given  for  writing  analysis 
codes.  The  data  organization  is  outlined  and  sample  data  is 
presented. 

This  publication  was  written  to  serve  as  an  addendum  to 
the  COPES  Manual;  "COPES — A  Fortran  Program  for  Engineering 
Synthesis,"  L.  E.  Madsen  and  G.  N.  Vanderplaats ,  NPS69-81- 
003,  Naval  Postgraduate  School,  Monterey,  California,  March 
1982. 
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II.  DESIGN  EXAMPLE 


It  is  required  to  design  the  cantilevered  beam  shown  in 
Figure  1.  The  objective  is  to  find  the  minimum  mean  volume 
of  material  which  will  support  the  concentrated  load  and 
maintain  an  allowed  probability  of  failure  of  10%. 


P 


i 


nA 

_ .  □  " 

L  =  200 !  N  — E>(  “£>j  B  |<I^ 


Figure  1.  CANTILEVERED  BEAM 


That  is , 


Minimize  volume  =  B  *  H  *  L  (1) 

The  mean  bending  stress  in  the  beam  must  not  exceed  its 
limit  which  is  20,000  psi  which  has  a  coefficient  of  variation 
of  0.07; 


BSTRES 


MC 

I 


20,000 


(2) 


The  mean  shear  stress  must  not  exceed  its  limit  of  10,000 
psi  which  has  a  coefficient  of  variation  of  0.05; 


SHRSTR  =  2X  =  755  -  10'000  (3) 

and  the  deflection  under  the  load  must  not  exceed  one  inch, 
coefficient  of  variation  equals  0.02; 

pr 3  4PL^ 

DELTA  *  £==■  =  <1.0  ( 4 ) 

3EI  EBH3  ~ 

Additionally,  geometric  limits  are  imposed  on  the  mean 
dimensions  so  that; 


0.5  1  B  <  5.0  (5) 

1.0  <  H  <  20.0  (6) 

I  <  10.0  (7) 

The  manufacturing  procedure  is  such  that  all  dimensions  have 
a  coefficient  of  variation  of  0.01.  The  mean  value  of  the 
dead  load  is  10,000  lb.  with  a  coefficient  of  variation  of 
0.07.  The  Young's  modulus  mean  value  is  30.E+06  psi  with  a 
0.06  coefficient  of  variation. 

The  ANALIZ  subroutine  on  page  11  of  the  COPES  manual 
is  used  for  the  analysis.  Only  one  modification  is  needed. 
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The  variable  P,  for  the  load,  is  added  to  the  global  common 

block. 

The  COPES  data  used  for  the  standard  optimization,  pages 
60  and  61  of  the  COPES  manual  is  modified  to  perform  proba¬ 
bilistic  optimization.  NCALC  in  data  block  B  is  changed  to 
7  to  indicate  probabilistic  optimization  is  to  be  performed 
and  IPROB  equals  1  is  entered  in  column  eight  so  the  normal 
distribution  model  will  be  used.  Three  additional  data 
blocks  are  required,  Ul,  U2  and  U3. 

DATA  BLOCK  Ul : 

There  are  five  probability  variables  which  make  up  the 
constraint  equations  2,  3,  4  and  7.  They  are  B,  H,  E,  AL, 
and  P;  so  IVAR  equals  5.  It  is  desired  that  gradients  of 
the  constrained  variables  only  be  calculated  at  the  beginning 
of  each  CONMIN  iteration,  therefore  LGRAD  equals  1  is  entered 
in  column  two. 

$  DATA  BLOCK  Ul 
5,1 

DATA  BLOCK  U2 : 

The  global  location  of  the  probability  variables  and 
their  corresponding  coefficients  of  variation  are  entered 
in  this  block. 

$  DATA  BLOCK  U2 


.01 

width,  B 

.01 

height,  H 

10, .05 


load,  P 

Young's  mod,  E 
length ,  AL 


8,  .06 
9,  .01 

DATA  BLOCK  U3: 

This  block  contains  the  allowed  probability  of  failure 
and  coefficient  of  variation  at  both  the  upper  and  lower 
bound  of  each  constrained  variable.  All  four  constrained 
variables  have  upper  bounds,  but  are  unbounded  on  the  lower 
end. 


$  DATA  BLOCK  U3 

1.  ,.07, . 1 , . 0  7 

BSTRES 

1. , .05, .1, .05 

SHRSTR 

1.  ,  .02, .1, .02 

DELTA 

1.  ,  .01,  .1, .01 

H/B 

The  allowed  probability  of  failure  at  the  upper  bound  and 
the  coefficient  of  variation  are  listed  in  columns  three 
and  four.  Because  there  are  no  lower  bounds  the  resulting 
probability  of  failure  will  be  zero.  These  first  two  columns 
will  be  ignored  during  prooability  of  failure  calculations 
in  COPES. 

Coupling  the  ANALIZ  subroutine  to  COPES  yields  the 
following  optimum  design: 
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CANTILEVERED  BEAM 


AL 

200.0 

P 

= 

0 . 10000E+0  5 

E 

= 

0 . 30000E+0  8 

B 

= 

1.92 

H 

= 

18.75 

VOL 

= 

7261.54 

BSTRES 

= 

0. 17802E+05 

SHRSTR 

= 

0.41727E+03 

DELTA 

= 

0 . 84  3  82E+00 

H/B 

= 

9 . 78 

The  COPES  output  for  this  example  is  provided  in  Figure  2. 


Ill .  PROGRAMMING  GUIDELINES 

The  programming  guidelines  presented  in  Chapter  II  of 
the  COPES  manual  all  apply  in  performing  probabilistic 
optimization.  One  additional  rule  must  be  followed  when 
writing  the  analysis  subroutine:  All  variables  in  the 
constraint  equations  which  have  probabilistic  distributions 
must  be  included  in  the  global  common  block  of  the  analysis 
subroutine. 


IV.  DATA  BLOCKS 


Input  instructions  for  the  probabilistic  optimization 


data  blocks  are  presented  in  the  following  five  pages. 
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V.  NOTES  TO  THE  USER 


1.  The  probabilistic  optimization  feature  is  not 
compatible  with  the  following  COPES  options: 

A.  sensitivity  analysis 

B.  two  variable  function  space 

C.  optimum  sensitivity 

D.  approximate  optimization 

2.  It  is  not  recommended  that  one  attempt  to  design  for 
allowed  probabilities  of  failure  of  less  than  0.00001  or 
greater  than  .99999  due  to  modifications  made  to  "widen"  the 
cumulative  normal  distribution  curve. 

3.  The  optimization  process  is  very  sensitive  to  the 
coefficients  of  variation  of  the  probability  variables  which 
make  up  the  constraint  equations.  Coefficients  of  variation 
of  less  than  0.0001  may  result  in  erroneous  results  if  the 
coefficients  of  variation  of  some  of  the  constraint  limits 
are  also  less  than  0.0001. 

4.  In  some  cases  the  initial  design  will  be  such  that 
the  optimization  process  terminates  unsatisfactorily  due  to 
the  calculation  of  zero  constraint  gradients.  Termination 
will  occur  either  due  to  an  inability  to  achieve  a  feasible 
design  or  CONMIN  will  reach  the  maximum  number  of  iterations 
no  matter  how  large  ITMAX  is  made.  This  probably  will  be 
due  to  the  fact  that  the  probabilities  for  failure  calcu¬ 
lated  are  either  1.0  or  0.0  and  the  finite  difference  steps 
are  insufficient  to  reach  the  region  along  the  cumulative 
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normal  distribution  curve  between  Pf  =  1.0  and  Pf  =  0.0. 

This  could  be  caused  by  an  initial  design  which  is  a  great 
distance  from  the  optimum  or  by  very  small  coefficients  of 
variation.  If  this  happens,  perform  a  deterministic  optimi¬ 
zation  to  achieve  a  design  in  the  vicinity  in  the  design 
space  desired  and  use  this  result  for  the  initial  design  in 
probabilistic  optimization. 

5.  In  order  to  determine  probabilities  of  failure,  COPES 
calculates  the  gradients  of  the  constrained  variables  with 
respect  to  the  probability  variables  which  make  up  the 
constraint  equations.  This  is  done  by  finite  difference 
steps  resulting  in  numerous  calls  to  the  ANALIZ  subroutine. 

In  many  problems  these  gradients  remain  essentially  the  same 
so  the  user  is  given  the  ability  to  determine  the  frequency 
in  which  these  gradients  are  calculated.  This  provides  for 
a  significant  reduction  in  the  number  of  function  evaluations. 
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